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NOMENC LA T URE 


% 


Symbol 

Definition 

Units 

A 

vent orifice area 

. 2 
in. 

CONF 

configuration number, identifies vent orifice 

— 

CORK 

correlate number, identifies one one 

plate POS, one CONF, and either one ratio 
Pp/Poo ° r ° ne P ro1:>e height h 

... 

h 

traversing probe height from flat plate 

in. 

K 

orifice efficiency (discharge coefficient) 

— 

M 

Mach number 

— 

P 

pressure 

psfa 

POS 

flat plate distance from tunnel wall 

in. 

q 

dynamic pressure 

psfa 

R 

gas constant, 1716 lb-sec/slug-°R 

— 

Re 

Reynolds number per unit length 

ft' 1 

RUN 

run number, identifies one M , one plate 
POS, and one CONF 00 

— 

SGFM 

standard volumetric flow rate 

ft^/min 

T 

vent orifice thickness 

in. 

T 

temperature 

°R 

V 

velocity 

ft/sec 

Greek 



y 

specific heat ratio, 1.4 

— 


viscosity 

lb -sec /ft' 

p 

density 

slug /ft ^ 
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Symbol 

Subscripts 

act 

DIT 

1 

j 

STD 
t, T 
th 
tr 
XL 
00 
1 

2 
9 


Definition 


actual (measured) 
quantity based on P 

oo 

index referring to pressures to P^q^ 
index referring to orifice jet 
standard conditions 
total conditions 
theoretical 

traversing probe or vicinity 
orifice extension lip 

r6fcrs to (or is based on) tunnel freestream 
conditions upstream of a shock wave 
conditions downstream of a shock wave 
refers to (or is based on) pressure P^ 
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Section 1 
INTRODUCTION 


During descent of space shuttle vehicles from high altitudes the build- 
up of ambient atmospheric pressures will result in crushing loads exerted 
across the vehicle skin. To relieve these loads and thereby circumvent struc- 
tural failure the vehicle compartment must be repressurized during descent 
such that the internal compartment pressures effectively counteract the exter- 
nal loads. It is anticipated that this equalization of pressures may be accom- 
plished by allowing external air to flow into the compartments through 
strategically located vents. Such vents have been used effectively in the 
Saturn/Apollo program to relieve potential bursting loads built up during 
ascent, where equalization of pressures was effected by allowing the trapped 
internal compartment gases to flow overboard to the reduced ambient pressure 
environment. The space shuttle vehicle, however, will require both compart - 
mental depressurization during ascent and repressurization during descent. It is 
currently anticipated that a vent-orifice system may be utilized during both 
phases of flight (Ref. 1). 

The prediction of venting performance of orifices under various external 
flow conditions requires a knowledge of orifice efficiencies for the conditions 
under consideration. Adequate outflow orifice efficiency data have been 
generated for flow conditions corresponding to the anticipated space shuttle 
ascent trajectory. These data were generated for application to the Saturn 
vehicle during a test program at NASA-Ames Research Center (Ref. 2). 

Existing inflow orifice efficiency data were limited, however, to a narrow range 
of external flow conditions. The test program described here was conducted to 
provide sufficient inflow data to enable venting analyses to be made for the 
space shuttle during the descent phase of flight. The pretest report associated 
with this test is presented in Ref. 3. 



Section 2 

TEST DESCRIPTION 



2.1 DESCRIPTION OF TEST FACILITY 

The test program was conducted at the NASA -Ames six-by- six-foot 
Supersonic Wind Tunnel, which operates on the closed circuit single-return 
principle. The throat area and location are controlled by an asymmetric, 
sliding-block nozzle. The test section of the tunnel is 14.4 ft long and 6 ft 
square in cross section. Tests can be performed at transonic conditions 
because of the perforated floor and ceiling of the test section. 

The test facility data acquisition system allowed raw data to be monitored 
during the test. The facility has eight vacuum pumps, two of which were used 
as the vacuum source to pull air through the vent plate orifices. The model 
could be viewed through two 46-in. diameter clear glass windows located on 
one side of the test section. The downstream window on the other side was 
removed and replaced by the flat plate model. 

2.2 MODEL INSTALLATION AND DESCRIPTION 

The model used for this test consisted basically of a flat plate. Fig. 1, 
which was traversed into the thick, turbulent boundary layer buildup on the 
tunnel wall. Maximum boundary layer thickness was obtained when the flat 
plate was flush with the tunnel wall. Minimum boundary layer thickness was 
obtained when the flat plate was traversed to a maximum position into the free- 
stream. Whenever the flat plate was in a position other than flush with the 
wall a sharp leading edge and side extensions attached to the plate minimized 
shocks and other undesirable flow effects created by the flow over the flat plate. 
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F ree stream 



Fig. 1 - Schematic of Model Setup 



The flat plate was instrumented with static pressure ports in lines 
parallel with the freestream and a total head rake positioned in line with the 
orifice. Figure 2 illustrates the static pressure port measurements obtained 
from (he flat plate. 

The rake probe measurements are shown numbered on Fig. 3. The rake 
was located 18.9 in. from the center of the vent orifice. A traversing probe 
was used to determine the boundary layer thickness at the location of the vent 
orifice. 

Five vent plates having various orifice geometries were used to com- 
pose the nine configurations shown in Fig. 4. The configurations were:a 0.75-in. 
diameter orifice with three depth configurations of 0.063, 0.15 and 0.3 in.; 
a four-to-one ratio ellipse with two depth configurations of 0.063 and 0.3 in., 
and one skew orientation configuration of 45 deg; one rectangular orifice with 
a length -to-width ratio of 18.12; one two-to-one ratio ellipse; and a 0.5 in. 
diameter circular orifice. The vent plates contained static pressure measure- 
ments on the plate surface and a lip static pressure port. The vent port 
designations are shown on Fig. 5. The lip pressure P is located just within 

the orifice, near ports P lr and P.-. 

r 15 40 

The vent plates were interchangeable with the plenum chamber. Instru- 
mentation in the plenum chamber included two static pressure measurements, 
one total pressure measurement and a temperature measurement. Downstream 
of the plenum chamber the flowmeter system, comprised of three flowmeters 
in parallel, recorded flow rate measurements. 

Reference 3 describes the test model and its measurements in detail. 

2.3 DATA ACQUISITION 

All data obtained during the inflow orifice efficiency test were initially 
recorded and reduced by the Ames facility. Subsequent analysis indicated that 
certain data, such as the bulk of flat plate and wake rake pressures, could be 
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Fig. 4 - Vent Orifice Descriptions 





Fig. 5 - Vent Plate and Port Pressure Measurement Designations 
(Idealized section shows addition of orifice extension.) 
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dropped from consideration as being largely insignificant and/or invariant. 
Similarly, certain other data, such as pressure parameters useful in pre- 
senting the orifice efficiencies, were added to the output. Thus the raw data 
were used to recompute a concise set of final data more immediately useful 
to venting studies than the original lengthy results. In addition, the actual 
mass flow rates were recomputed using a more accurate table of instrument 
calibration data. 

Notable changes from the scheduled output of the pretest report (Ref. 3) 
are as follows: 


• Plenum pressure is always P . Only minor 

P Pi 

differences were found between P and P . 

P 1 p 2 

• Plenum total pressure P fc was dropped due to the 

p 

continuously varying flow direction within the plenum 
chamber . 


oo 


• The theoretical plenum (orifice) temperatures T fch 

G 

and T were not directly considered during the 
th 9 

analysis, and therefore dropped. 

• A third orifice efficiency, K DIT * was added to the 
final results. This uses the assumptions of Ref. 4. 


• Flat plate, wake rake, and traversing probe pressure 
data are represented in the figures of Appendixes A, 
B, and C. They are deleted from the data listing. 


2.4 DATA REDUCTION EQUATIONS 

Most basic data are obtained through direct conversion of the various 
instrument outputs. The derived data are obtained from the equations de- 
scribed on the following pages. 
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2.4.1 Boundary Layer Traversing Probe and Wake Rake Pressures 


The total pressure readings for both the probe and rake are corrected 

for standing shock waves when the freestream flow is locally supersonic. 

With subscripts 1 and 2 representing conditions upstream and downstream of 

the shock, respectively, the criterion for this correction is P./P. < .5283. 

% * *2 
Otherwise no correction is needed. 

* The solution involves solving the Rayleigh pitot formula, 



(y + 1) Mj 

y-1 

Y + 1 

P ! ' 

2 

vi 


_2y Mj - (y-1) _ 


1 

y-1 


for Mj ahead of the shock; then 


P 


t 


1 




" (y + 1) Mj 

y-1 

y+1 

(y-1) Mj + 2 


_2y Mj - (y-l)_ 


( 1 ) 


( 2 ) 


An approximate explicit solution of Eq. (1) was used in the test data reduction. 

For the traversing probe, the following substitutions are made in Eqs. (1) 
and (2): 

Pj = = local upstream static pressure 

P t = (P tr )2 = probe reading 
2 

Pj. = P^ r = corrected pressure. 
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For the wake rake, these substitutions are 


P = P_„ (located 3.9 in. upstream of rake on the flat plate) 

local upstream static pressure 

P = (P.) ? = i** 1 rake tube reading 
t 2 1 L 

P = P. = corrected pressure 
t l 1 


2.4.2 Traversing Probe Local Velocity 


V 


tr 



ft/sec 


where P. = corrected pressure from above, 
tr 

2.4.3 Local Flat Plate Mach Number (based on P^) 


hJ. 



2.4.4 Actual Mass Flow Rate 


(pAV) act = < SCFM > slug/sec 


where SCFM = corrected flowmeter output interpolated on calibration 
curves, ft^/min 

pSTD = .0023275 slug/ft 3 


2.4.5 Mass Flow Ratios 


These are ratios of mass flow rate per unit area, where the numerator 
is 




* 


(pv). . 

<pA v > act /A 

(6) 

of the two values 


- 

P V /RT 
00 00 00 

(7) 

p ,m 9 [ 

nf- <• ♦ ¥ ^ >r 

(8) 

L 

00 J 



2.4.6 Orifice Efficiencies 

The orifice efficiency is obtained as the ratio of actual (measured) 
to theoretical (ideal) mass flow rate: 


(p Av ) act 

(pAV) th 


(9) 


Actual flow rate is computed in Eq. (5). Theoretical flow rates are 
determined by applying the relations for steady, one -dimensional, compres- 
sible flow of a perfect gas to the measured conditions across the vent orifice. 
Several choices are available for the total pressure P T through the orifice; 

in all cases the total temperature is set equal to the freestream total tempera 
ture L . 


00 
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Two limiting conditions are applied to Eq. (10): 

For choked flow, where P p /P T < 0.5283, the theoretical flow rate is 


<P AV >th 


A P. 



y + l 

(0.5283) y 


R T, 


00 


(1 


When the ratio P /P T > 1.0, 

P T — 

(pAV) th = 0.0 < 

and the orifice efficiency K is also set equal to zero. 

These efficiencies were computed with the following choices for total 
pressure P^: 

• P^ = Pg, and K = Kg. This choice follows Ref. 5. P 9 is a flat 

plate static pressure located 7.5 inches upstream of the 
orifice. 

. P T . P„. and K = K,. 

• p_ = P. , and K = K___. This follows Ref. 4. 



Section 3 

PRESENTATION AND DISCUSSION OF DATA 

This section contains a description of the data shown in Appendixes A 
through E and a discussion of the results. The variables plotted are representa- 
tive of the overall objectives of the test requirements. Several modifications 
were made to the test program as presented in the pretest report (Ref. 3). 

The major modification was the deletion of two vent plate configurations (6 
and 7) and the inclusion of vent plate configuration 11. The latter modification 
was made during the course of the test. The orifice area of configuration 11 
was sufficiently small (0.1964 in. 2 ) for the Ames pump system to attain the 
desired pressure drop across the orifice. 

3.1 DATA PRESENTATION 

The test data are presented in Appendixes A through E in plotted form 
and in Appendix G in tabulated form. Appendix F presents the nomenclature 
used to define the data in Appendix G. 

3.1.1 Traversing Probe 

The traversing probe was used to measure total pressure at various 
distances from the flat plate. The total pressure measurements were con- 
verted to velocities using the Rayleigh pitot tube formula. Figures A-l 
through A -4 in Appendix A are results of the velocity distributions normalized 
to the freestream velocity for different Mach numbers compared at the three 
flat plate positions 0.0, 1.75 and 5.85 in. The boundary layer was assumed 
to be fully developed at 99% of the freestream velocity. The resulting boundary 
layer thickness as a function of Mach number for each plate position is pre- 
sented in Fig. A -5. 


3.1.2 Wake Rake 


Wake rake pressure distributions representative of the entire range of 
flows encountered during the test are presented in Appendix B in Figs. B-l 
through B-6. The rake, located 18.9 in. ait of the vent area had a total of 23 
pitot tubes at heights ranging from 0.04 to 7.0 in. above the flat plate. Each 
rake tube pressure was nondimensionalized by the freestream total pressure. 

3.1.3 Plate Static Pressure Distributions 

The plots in Appendix C, Figs. C-l through C-5, describe the static 
pressure distributions on the longitudinal axis of the flat plate model in line 
with the vent orifice. The static pressures are normalized by the tunnel 
freestream total pressure. The pressure distributions in Appendix C are 
representative of the maximum suction through the orifices of configurations 
1, 2, 9 and 1 1 at the maximum flat plate boundary layer thickness; i.e., POS = 
0.0 in. The pressure distribution associated with configuration 9 on Fig. C-5 
is also plotted at the minimum flat plate boundary layer thickness; i.e., POS 
5.85 in. 

3.1.4 Orifice Efficiencies 

Figures D-l through D-54 presented in Appendix D are the results of 
the orifice efficiency test data. The curves on Figs. D-l through D-27 show 
the orifice efficiency K^, based on local pressure P^, versus the local to 
plenum static pressure ratio P^/P^. Each plot includes all Mach numbers 
run for the specific configuration and flat plate position. The next set of 
curves. Figs. D-28 through D-54, shows the jet to local mass flow ratio, 
(pV)./(pV)g), versus the local to plenum pressure ratio for these same con- 
figurations . 


3.1.5 Configuration Comparisons 


Other requirements of the in-venting efficiency test were to determine 
the effects of varying the vent orifice geometry, vent thickness, vent 
orientation, and boundary layer thickness. Figures E-l through E-20 of 
Appendix E are results of such comparisons. Figures E-l through E-3 show 
the effects of variation between the size of circular orifices (0.5 and 0.75 in. 
diameter) at Mach numbers of 0.7, 1.1 and 1.9 at plate positions of 0.0, 1.75 
and 5.85 in. Vent plate thickness effects at various Mach numbers are shown 
on Figs. E-4 through E-6 for the circular orifice and on Figs. E-7 through E-9 
for the orifice of aspect ratio 4:1, Vent orientation with respect to freestream 
flow direction is compared in Figs. E-10 through E-12 for the vent orifice of 
aspect ratio 4:1. The next configuration comparisons are aspect ratio effects 
in Figs. E-l 3 through E-l 5. Compared are aspect ratios of 1 : 1, 2:1 and 4:1 
for different Mach numbers at plate positions 0.0, 1.75 and 5.85 in. Finally, 

Figs. E- 16 through E-20 show effects of plate position for several configurations. 

3.2 DATA DISCUSSION 

The following is a discussion of the overall data trends noted in the test 
results. 

3.2.1 Traversing Probe Data 

Data obtained from the traversing probe tests provided the information 
to define boundary layer thickness over the flat plate as a function of free- 
stream Mach number for the plate positions investigated during the test. 
Inspection of Figs. A-l through A-4 reveals a consistent variation of boundary 
layer thickness as a function of plate position and Mach number. The results, 
plotted in Fig. A -5, are compared with the out-flow venting results of Walters 
(Ref. 2). Good agreement is noted in the comparisons for positions 0.0 and 
5.85 in. The wavy trend in the outflow data is not indicated in the inflow 

data; this fact may be due to the reduced quantity of Mach number data points 
in the inflow data. 



3.2.2 Wake Rake Data 


Wake rake measurements of total pressures were included in this test 
program to provide some indication of the effect of inflow venting on the 
boundary layer profile downstream of the orifice. Varying the orifice inflow 
resulted in negligible perturbation of the boundary layer, indicating complete re- 
establishment of external flow at the rake. Representative data from Configura- 
tion 1 are shown in Figs. B-l through B-6 for various Mach numbers and plate 
positions. 

3.2.3 Plate Static Pressure Data 

The effect of inflow venting on static pressure distribution on the flat 
plate is presented in Appendix C. The dip in the trend prior to the vent port 
indicates a velocity increase, while the peak following the vent port indicates 
a velocity decrease. Both the accelerated flow upstream of the orifice and 
the decelerated flow downstream are due to suction through the orifice. Further 
downstream the flow returns to the freestream velocity. The static pressure 
is seen to increase beginning at about 15 in. downstream of the orifice for 
Mach numbers 1.1 and 1.3. This is probably due to the shock located at the 
boundary layer rake 18.9 in. downstream of the orifice. 

3.2.4 Orifice Efficiency Data 

The orifice efficiency data are presented in two forms in Appendix D. 

The first form shows orifice efficiency plotted versus the ratio of local pres- 
sure to plenum pressure using Mach number, configuration, and plate position 
as parameters. The orifice efficiency is based on local pressure P^, located 
7.5 in. upstream of the orifice. 

The second set of orifice efficiency data was plotted in the form of jet- 
to-local mass flow ratio versus local-to-plenum pressure ratio. These curves 
were included because of their overall smoother trend as compared to the 
orifice efficiency K^. 
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It was noted that the data curves do not generally achieve zero efficiency 
or mass flow ratio at a pressure ratio of 1.0. This was especially evident in the 
data based on P. as compared to that based on indicating a local boundary 
layer effect. For this reason, no curves of are included here. 

3.2.5 Configuration Comparison Data 

Appendix E shows the results of comparing configurations for the effect 
on orifice efficiency due to varying orifice size, vent plate thickness, vent plate 
orientation, orifice aspect ratio, and plate position. The comparisons were 
made with mass flow ratio for Mach numbers 0.7, 1.1 and 1.9 at plate positions 
of 0.0, 1.75 and 5.85 in. from the tunnel wall. 

A comparison in Figs. E-l through E-3 of the results for circular orifices 

(0. 5 and 0.75 in. diameter) shows that the larger orifice provides the greater ef- 
ficiency. 

A comparison of the circular vent plate thicknesses presented in Figs. E-4 
through E-6 shows, in general, that the thinner vent plate is the most efficient. 
This trend is consistent with one exception on Fig. E-4 at Mach 1.9 and with the 
maximum boundary layer thickness (position 0.0 in.). The anticipated result is 
est expressed by Fig. E-5 in which the thickest of the vent plates is the least 
efficient, while the thinnest vent plate is the most efficient. This same trend is 
also noted in the elliptical vent plate comparisons of Figs. E-7 through E-9. The 
result is even more obvious as the Mach number increases. 

The comparisons on Figs. E-10 through E-12 are based on orientation of 
the elliptical vent plate (aspect ratio of 4:1) with respect to the freestream air 
flow. Orientation of the major axis of the ellipse parallel to the flow direction 
resulted in higher efficiencies than the 45 deg orientation. No data were taken 
at intermediate orientation positions. The difference in efficiency for the two 
orientations increased with Mach number. 

The next set of comparisons (Figs. E-13 through E-15) was based on the 
effects of orifice aspect ratio on inflow. The curves obtained from orifices of 
approximately equal area and aspect ratios of 1:1, 2:1, and 4:1 show a substantial 
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increase in efficiency with increasing aspect ratio (major axis parallel to flow 
direction) and decreasing Mach number. 

Figures E-16 through E-20 show effects due to varying flat plate posi- 
tion on configurations 11, 1, 3, 2, and 9. It is seen that varying the plate 
position has a definite effect on efficiency, especially for the smaller orifices. 
The lack of consistent trends, however, indicates further analysis is warranted. 
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Appendix A 

TRAVERSING PROBE DATA 
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Fig. A-l - Velocity Ratio Profile at Mach 1.891 
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Fig.A-4 - Velocity Ratio Profile at Mach 0.700 
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Fig. A- 5 - Comparison of Boundary Layer Thickness Over the Test Plate at the Ames 

6x6 Foot Wind Tunnel as Obtained from the Inflow Venting Test and from Ref. 
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Appendix C 

PLATE STATIC PRESSURE DATA 
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Fig. C-3 - Static Pressure Effects of Inflow Venting at Different Mach Numbers 


Position 0.0 in. Configuration 9 
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Fig. C-4 - Static Pressure Effects of Inflow Venting at Different Mach Numbers 
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Fig.C-5 - Static Pressure Effects of Inflow Venting at Different Mach Numbers 
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ORIFICE EFFICIENCY DATA 
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Orifice Efficiency vs Pressure Ratio for Configuration 11 and 
Plate Position 0.0 Inches for Various Mach Numbers 
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D- 10 


I 




D- 1 1 


o no m 


o r^ 
-I o o 



o 

ro 


00 


o< 

vO \ 

o 

^ CX» 


^ (TJ 

oS 

0) 

1 5 

(W 

(A 

a> 

on ^ 

_; CU 


(M 


O' 


D- 12 


Fie D-12 - Orifice Efficiency vs Pressure Ratio for Configuration 10 and 
Plate Position 5.85 Inches for Various Mach Numbers 
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Fig. D- 16 - Orifice Efficiency vs Pressure Ratio for Configuration 2 and 
Plate Position 0.0 Inches for Various Mach Numbers 
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Orifice Efficiency vs Pressure Ratio for Configuration 2 and 
Plate Position 5.85 Inches for Various Mach Numbers 
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Fig. D-19 - Orifice Efficiency vs Pressure Ratio for Configuration 5 and 
Plate Position 0.0 Inches for Various Mach Numbers 
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Fig. D-23 - Orifice Efficiency vs Pressure Ratio for Configuration 8 and 
Plate Position 1.75 Inches for Various Mach Numbers 



Fig. D-24 - Orifice Efficiency vs Pressure Ratio for Configuration 8 
and Plate Position 5.85 Inches for Various Mach Numbers 
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Fig. D-25 - Orifice Efficiency vs Pressure Ratio for Configuration 9 
and Plate Position 0.0 Inches for Various Mach Numbers 
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Fig. D-26 - Orifice Efficiency vs Pressure Ratio for Configuration 9 
and Plate Position 1.75 Inches for Various Mach Numbers 
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Fig. D-27 - Orifice Efficiency vs Pressure Ratio for Configuration 9 
and Plate Position 5.85 Inches for Various Mach Numbers 
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Fig. D-42 - Mass Flow Ratio vs Pressure Ratio for Configuration 3 
and Plate Position 5.85 Inches for Various Mach Numbers 
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Fig. D-43 - Mass Flow Ratio vs Pressure Ratio for Configuration 2 
and Plate Position 0.0 Inches for Various Mach Numbers 
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Fig. D-44 - Mass Flow Ratio vs Pressure Ratio for Configuration 2 
and Plate Position 1.75 Inches for Various Mach Numbers 
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Fig. D-49 - Mass Flow Ratio vs Pressure Ratio for Configuration 8 
and Plate Position 0.0 Inches for Various Mach Numbers 
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Fig. D-51 - Mass Flow Ratio vs Pressure Ratio for Configuration 8 
and Plate Position 5.85 Inches for Various Mach Number 
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Fig. D-54 - Mass Flow Ratio vs Pre 
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CONFIGURATION COMPARISON DATA 
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Fig. E-2 - Comparison of Mass Flow Ratio of Circular Orifice 

Configurations with Different Areas at Plate Position 
1.75 Inches 
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Fig. E-4 - Comparison of Mass Flow Ratio of Circular Orifice 
Configurations with Different Thickness at Plate 
Position 0.0 Inches. 
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Fig. E-7 - Comparison of Mass Flow Ratio of Elliptic Orifices 

of Aspect Ratio 4:1 with Different ThicKness at 
Plate Position 0.0 Inches. 


E - ' 














r 


POS 5.85 in. 


Mach 


CONF2 


1.9 


0 


1.1 


0.7 


A 


CONF8 












Mass Flow Ratio (pV)j /(pV)g 



Fig. E-10 - Comparison of Mass Flow Ratio of Elliptic Orifices 

at Different Orientations at Plate Position 0.0 Inches. 
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Fig. E-15 - Comparison of Mass Flow Ratio of Elliptic Orifices 

with Different A spect Ratios at Plate Position 5.85 
Inches. 
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Fig. E-17 - Comparison of Mass Flow Ratio for Configuration 1 
at Different Plate Positions 
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Fig. E- 19 - Comparison of Mass Flow Ratio for Configuration 2 
at Different Plate Positions 
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Appendix F 
OUTPUT NOTATION 



Symbol 

Definition 

Units 

CORR 

correlate number; identifies one POS, one 
CONF, one M, one PP/P 



CONF 

configuration number; identifies vent plate 

— 

POS 

flat plate position from tunnel wall 

in. 

AREA 

area of vent plate orifice 

. 2 
m 

PT 

freestream total pressure 

lb/ft 2 

P 

freestream static pressure 

lb/ft 2 

TT 

freestream total temperature 

°R 

T 

freestream static temperature 

°R 

M 

freestream Mach number 

— 

Q 

freestream dynamic pressure 

lb/ft 2 

V 

freestream velocity 

ft/sec 

RU 

freestream Reynolds number per foot 

ft -1 

P9 

flat plate static pressure Pq, located 7.5 
inches upstream of orifice ' 

lb/ft 2 

P9/P 

local to freestream static pressure ratio 

— 

M9 

local Mach number based on P9 

— 

P16 

orifice lip static pressure 

lb/ft 2 

PXL 

orifice extension lip static pressure 

lb/ft 2 

TP 

plenum chamber static temperature 

°R 

tp/tt 

plenum static to freestream total temper- 
ature ratio 


PP 

plenum chamber static pressure 

lb/ft 2 

MDOT 

mass flowrate through orifice 

slug/ sec 

MVJ 

orifice unit area mass flowrate 

slug/ft 2 - 

MVINF 

freestream unit area mass flowrate 

slug/ft 2 - 

MV9 

local unit area mass flowrate 

slug/ft 2 - 

KINF 

inflow orifice coefficient based on free- 
stream static pressure 

— - 

K9 

coefficient based on local static pressure 

— 

KDIT 

coefficient based on freestream total pressure 

— 


F-l 


sec 

sec 

sec 
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NOMENCLATURE 


Svmbol 


Units 

A 

vent area 

. 2 
in 

C 

p 

specific heat at constant pressure 

Btu/lbm-mole-°R 

C 

V 

specific heat at constant volume 

Btu/lbm-mole-°R 

n 

enthalpy of gas 

/l — w rtl £* 

J UTVlAf 

K 

orifice efficiency 

dimensionless 

M 

Mach number 

dimensionless 

m 

mass of gas 

slug 

m 

mass flow rate 

slug/ sec 

P 

pressure 

lb/ft 2 

P 

time rate of change of pressure 

lb/ft 2 /sec 

R 

gas constant (1716 for air) 

ft-lb/ slug-°R 

t 

time 

sec 

T 

temperature 

°R 

V 

compartment volume 

ft 3 

Greek 



Y 

P 

ratio of specific heats = C^/C^ 
gas density 

dimensionless 

slug/ft 3 

Subscripts 



c 

» 

denotes compartment 


P 

denotes compartment or plenum chamber 


1 

denotes starting condition of time step 


2 

denotes end condition of time step 


o 

denotes total conditions 


00 

denotes local external conditions 


j 

denotes orifice jet Conditions 




Section 1 
INTRODUCTION 


During descent of space shuttle vehicles from high altitudes the buildup 
of ambient atmospheric pressures will result in crushing loads exerted across 
the vehicle skin. To relieve these loads and thereby circumvent structural 
failure, the vehicle compartments must be repressurized during descent such 
that the internal compartment pressures effectively counteract the external 
loads. It is anticipated that this equalization of pressures may be accomplished 
by allowing external air to flow into the compartments through strategically lo- 
cated vents. Such vents have been used effectively in the Saturn/Apollo program 
to relieve potential bursting loads built up during ascent, where equalization of 
pressures was effected by allowing the trapped internal compartment gases to 
flow overboard to the reduced ambient pressure environment. The space shut- 
tle vehicle, however, will require both compartment depressurization during 
ascent and repressurization during descent. It is currently anticipated that 
a vent -orifice system may be utilized during both phases of flight (Ref. 1). 

The prediction of venting performance of orifices under various external 
flow conditions requires a knowledge of orifice efficiencies for the conditions 
under consideration. Adequate outflow orifice efficiency data have been gener- 
ated for flow conditions corresponding to the anticipated space shuttle ascent 
trajectory. These data were generated for application to the Saturn vehicle 
during a test program at NASA-Ames Research Center (Ref. 2). Inflow orifice 
efficiency data have been limited, however, to a narrow range of external flow 
conditions. 

The purpose of the investigation described in this report is to provide 
sufficient inflow orifice efficiency data and to develop the computer program 
capability to enable compartment venting analyses to be made for the space 
shuttle durine the descent phase of flight. 



Section 2 

INFLOW ORIFICE EFFICIENCY TEST 
2.1 TEST DESCRIPTION 

A test program was undertaken to determine orifice efficiencies for the 
flow of air into a compartment from a flowing external stream. The test pro- 
gram was conducted at the NASA-Ames 6 x 6-Foot Supersonic Wind Tunnel 
during a two-week period (20 September 1971 through 1 October 1971). This 
particular wind tunnel was selected because hardware was available from a 
previous test program for outflow venting which was conducted in this same 
wind tunnel (Ref. 2). The hardware was designed to fit into the window of the 
Ames 6 x 6-foot tunnel, and only minor modifications of the test hardware were 
required. Also, the Mach number range capability of the tunnel, 0.6 to 2.2, 
fell within a major portion of the anticipated shuttle descent trajectory Mach 
number range (Ref. 1). 

The test model consisted basically of a flat plate which was designed to 
traverse into a boundary layer created by the turbulent air buildup on the tunnel 
wall. Other pertinent components of the model include a vent plenum, a 
traversing probe plenum and interchangeable vent plates with variously shaped 
orifices. The vent plenum with the interchangeable vent plates allowed for 
air removal by means of a vacuum system. The measured mass flow through 
the orifice and pressure measurements on the flat plate and plenum chamber 
supplied information necessary to calculate the orifice flow coefficient. The 
traversing probe plenum supplied information on boundary layer thickness over 
the vent port location. 

A Mach number range of 0.7 to 1.9 was used during the test along with a 
tunnel total pressure of 1584 psfa. Other parameters varied were orifice 
geometry, vent plate thickness, vent orientation, flat plate boundary layer 
thickness and pressure ratio across the vent plate. 
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Data obtained from the test include three calculated inflow orifice effi- 
ciencies based on freestream static and total pressure and a plate static pres- 
sure port located 7.5 inches upstream of the vent port. Other data obtained 
include traversing probe data to define boundary layer thickness, boundary 
layer rake data at a location 18.9 inches downstream of the vent port, and flat 
plate static pressure data to define the longitudinal pressure distribution along 
the centerline of the flat plate. This test program and the experimental results 

* are described in detail in the pretest and test reports (Refs. 3 and 4). 

2.2 TEST RESULTS AND DISCUSSION 

Inflow orifice efficiencies were obtained as three distinct coefficients 
differing in the choice of total pressure through the orifice. These were: 

• K^, based on local static pressure Pg (7.5 inches upstream of 
orifice). This follows the method of Kalivretenos, Ref. 5. 

• K , based on freestream static pressure, . 

• based on freestream total pressure, P^. This 
follows the method of Dittrich, Ref. 6. 

A comparison of Kg and K # is shown in Fig. 1 for a circular orifice of 
0.75 inch diameter and a maximum boundary layer thickness (plate position 
0.0). The curves of constant Mach number are uniformly spaced and con- 
vergent for Kg, but shift relative to each other for K oQ . This was due to 
appreciable differences between the static pressures Pg and P^ in general. 
Figure 2 shows the ratio Pg/P ^ as a function of Mach number and plate posi- 

* tion for several orifice configurations, obtained for zero or near -zero inflow 
venting. These unexpected large differences are not easily accounted for, but 
possible contributing factors may include the plate supporting structure and 
blockage alleviation system being located on the same wall as, and downstream 
of, the pressure port used to determine P^. 

In any case, the consistently smoother trends for Kg and other param- 
eters based on Pg encouraged data analysis to be concentrated on results 
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using this local pressure rather than 

Results of this test are compared with those of the two major inflow 
venting references used, Kalivretenos and Dittrich, in Figs. 3 through 6. 
The notations of both are changed to the common form used here. Figure 3 
reproduces the data of Kalivretenos for various circular orifices. The co- 
efficients are plotted against a flow parameter, (P n -P )/(P, -P 0 ), where 

7 p t co 7 

P is the compartment static pressure. Results of the present study are 
H 

~ ~ j _ tt 1 z ~ a : ~ ~ i i- xi - r r\ j 

in x ig. -x x ii a. oiiuiiax laoiiiua, aiuuuu^u LHC I'diigt: ui uuw paiamcici 

achieved is much lower due to reduced absolute pressure levels and equip- 
ment limitations. 


Figure 5 reproduces the Dittrich data for a circular orifice where the 

coefficient is based on total freestream pressure, and the flow parameter 

is (P, -P ) /( P*. -Pq). Results of the present test are shown in Fig. 6; again, 

^ 00 P 00 ' 

the flow parameter range is small. 

These figures show some marked differences. In Figs. 3 and 4, for 
example, the trends with increasing Mach number for a constant value of 
flow parameter are exactly opposite. Figures 5 and 6 indicate a similar 
discrepancy. 

Differences are also noted in general conclusions. Both references 
conclude that boundary layer thickness effects on orifice efficiency are neg- 
ligible while this study found considerable variation with plate position, if 
not boundary layer thickness directly; however, the boundary layer in Dittrich 
was on the order of 0.1 inch thick while present results are for thicknesses 
up to 4.5 inches. 

No analysis of these discrepancies was attempted here. Known differ- 
ences among the three tests were total pressure magnitude, boundary layer 
thickness, sheer facility size, and possible instrument error. The incon- 
sistencies appear sufficiently great to warrant further study. 
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Orifice Orifice 
Diam,(in.) Thickness (in.) 
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Fig. 3 - Orifice Efficiency vs Flow Parameter for Circular Orifices, 
from Kalivretenos et al. (Ref. 5) 




Orifice Efficiency vs Flow Parameter for 0.75 in. Diameter Circular Orifi 
from Dittrich et al. (Ref. 6). Orifice Thickness is 0.040 in. 








Figure 7 shows another parameter chosen to represent venting effi- 
ciency, the ratio of orifice to local external mass flow rate, (pV)./(pV) 
This parameter was also used in presenting the test results in Ref. 3. 

A summary of the data trends from the inflow test program gives the 
following results in terms of flow efficiency: 


• Orifice efficiency is primarily a function of external flow velocity 
and pressure ratio across the orifice. It increases with decreas- 
ing external velocity and with increasing external -to -chamber 
pressure ratio. 

• Efficiency increases with increasing orifice area. 

• Efficiency increases with decreasing wall thickness. 

• Efficiency increases with increasing orifice aspect ratio (long 
axis aligned with external flow direction). 

• Efficiency is maximized for orifice aspect ratios greater than 
one when the major axis is aligned parallel with external flow 
direction. 

• Efficiency is affected by boundary layer thickness, but the 
degree of effect has not been analyzed. It is noted that in- 
creasing orifice size reduces the effects. 


The effects on orifice efficiency due to geometry are probably linked 
with viscous phenomena and blockage at the downstream orifice edge. An 
entering stream tube is constricted by the upstream boundary layer and the 
downstream orifice wall; only the latter boundary can be controlled. Geom- 
etry changes such as decreased wall thickness and increased flow entry length 
(by a larger, longer, or more closely aligned orifice) all contribute towards 
a larger entry area and hence more efficient flow. 


* 
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Section 3 


OUTFLOW VENTING COMPUTER PROGRAM 


A computer program originally developed by Jump and Henson (Ref. 7) 
serves as the basis for both the inflow and outflow venting programs used in 
this report. Basic features of the outflow program are described below, as 
are results for a typical space shuttle ascent trajectory. 

3.1 OUTFLOW COMPUTER PROGRAM DESCRIPTION 

The procedure used in Ref. 7 to calculate compartment pressure time 
histories from vehicle launch is a numerical operation. The program pro- 
duces calculations for given time steps beginning from launch until the atmo- 
spheric ambient pressure reaches essentially vacuum conditions. 

Compartment pressures are calculated by use of a numerical procedure 
which works as follows: 

• The change of mass in the compartment between time 
steps is calculated from the isentropic relation p c /p ( 7 = 

constant after "guessing” the compartment pressure at 
the end of the time step. 

• The calculated change of mass is balanced against the 
average mass flow rate through the vent port for the 
same time interval as established by the continuity 
equation, determining a second value of compartment 
pressure. Outflow orifice efficiencies obtained from 
the data of Ref. 2 are used in these computations. 

• If the pressures are equal, within an allowed error 
range, then the program advances to the next time 
step. 

• If the pressures are not equal, then a new compartment 
pressure is "guessed" and the procedure is repeated. 


3.2 OUTFLOW COMPUTER PROGRAM RESULTS 

A typical space shuttle ascent trajectory is illustrated in Fig. 8. The 

outflow program was run using this trajectory information and outflow orifice 

efficiency data for a 4:1 ratio elliptical orifice. Results of severalsuch runs 

are shown in Fig. 9. The graph shows the skin differential pressure as a 

3 . * 

function of time from launch for a compartment volume of 3000 ft and orifices 

2 

of area 50, 75 and 150 in . These calculations were based on the assumption 
that the local external pressure is identical to the freestream ambient pressure. % 
The curves show that an orifice area to compartment volume ratio of 0.0166 
in /ft must be maintained to keep the skin differential pressure within as- 
sumed design limits of 1.0 psi. (Ref. 1). 
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Fig. 8 - Typical Space Shuttle Booster Trajectory Data 




Skin Differential Pressure as a Function of Time for a Typical. Ascent Trajectory 


Section 4 

INFLOW VENTING COMPUTER PROGRAM 

4.1 INFLOW COMPUTER PROGRAM 

The basic outflow venting computer program written by Jump and 
Henson (Ref. 7) was modified to calculate inflight compartment pressures 
during inflow venting. The program was designed to handle one compart- 
ment with up to four vent orifices. The present version of the program 
does not allow for outflow venting, and hence is limited to descent trajec- 
tories. Future versions will incorporate both inflow and outflow capabilities. 

The program utilizes a numerical calculation procedure in the solution 
of the following equations: 

The equation of state applied to the gas inside the compartment is 

P = -§-m T (1) 

c V 

The time rate of change of the compartment pressure, therefore, is 

p c = v- -ar (mT) (2) 

The rate of energy flow into the compartment is the product of the mass flow 
rate and the local total enthalpy. This energy flow rate is equated to the time 
rate of increase in the compartment internal energy as follows: 

- h o = C V -$£- ( mT ) (3) 


Combining Eqs. (2) and (3) and employing the definition of total enthalpy 

as h = C T yields the rate of change of the compartment pressure: 
o p o 7 



( 4 ) 
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The numerical calculation corresponding to Eq. (4), assuming linear varia- 
tion of all parameters over a small time increment At, is 


P = P + At 
c 2 C 1 



’ T + T 


• * 

°1 °2 


m i + m ^ 

2 


2 


(5) 


where the subscripts 1 and 2 refer, respectively, 
and end of the time interval At. The values of T q 


4b a #1 *1 ^ < /*% O 

tu 


-> +■ +V-» Vwanri na 

«■*- . — J~> 


are those of the local external 


total temperature. 
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Parameters known at the beginning of the time step include the following: 


P c , At, V, R, V, T , T 0 , P„ 
1 “ 


m. 


'1 


A tentative calculation of m 2 is made by assuming a value for the compartment 

pressure P equal to P . 

c 2 1 


The mass flow through the orifice at the end of any time step At is 



( 6 ) 


The term K 2 is the experimentally obtained orifice efficiency at the end of 

At and is a function of P /P and M . A more detailed explanation and 

c 2 00 2 °°2 

derivation of Eq. (6) is presented in Refs. 3 and 4. 
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The assumed value of P is used in Eq. (6) to determine a flow rate 

r ^ l 2 w ^ich then substituted in Eq. (5) to recompute a new value P' . The 

c 2 

iteration cycle reduces the difference between F c and P' to an allowable 
error, and the next time step is begun. ^ ^ 

When there are several mass flow sources (vents), Eq. (6) is employed 
separately for each, the mass flow rates m ^ and of Eq. (5) are the sum - 
mations over all sources. The values of T 0j and T <> 2 in Eq. (5), however, 
are weighted averages of all sources according to orifice area. 

Appendix A is a user's manual describing the applications, restrictions, 
and input requirements of the program including a listing of sample input data. 
Appendix B is a listing of the complete inflow program. 
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4.2 INFLOW COMPUTER PROGRAM RESULTS 

A set of compartment pressure histories was generated using the in- 
flow computer program with inputs of the booster descent trajectory of Fig. 8 
and experimental orifice efficiency data. These computed pressure histories 
are presented in Fig. 10 for a compartment of 3000 ft 3 volume and for elliptical 
vents of 4:1 aspect ratio with areas of 10 and 50 in . 

As in the outflow venting case, the local external flow conditions are 
assumed to be identical to freestream conditions. The vent area to compart- 
ment volume ratio of 0.0166 in 2 /ft 3 (50 in 2 orifice) applied to the descent 
trajectory shows a peak pressure differential of only 0.2 psi as compared to 
the 1.0 psi for the same vent orifice in the outflow venting case. The orifice 
area of 10 in 2 is representative of a typical minimum equivalent leakage area 
for a spacecraft compartment volume of 3000 ft 3 . The minimum equivalent 
leakage area was obtained from Saturn flight vehicle AS- 504 data as published 
in Ref. 8. The interstage compartment equivalent leakage area was deter- 
mined by pressurizing the vehicle and monitoring the compartment pressure. 
The area corresponding to the equivalent leak rate was then termed the equiv- 
alent leakage area. As seen in Fig. 10, the maximum pressure differential of 
the compartment, assuming a minimum equivalent leakage area, is about 1 psi. 

Figure 11 traces the compartment pressure histories during descent for 
vent areas of 10 and 50 in 2 . Figure 12 shows the corresponding compartment 
temperature histories. Note that a peak compartment temperature of about 
350°F occurs at about 520 sec descent flight time for the 50 in orifice. After- 
ward, the temperature drops to a level of about 150 to 17 5°F which it maintains 
during the remaining descent. At the time of peak compartment temperature, 
the compartment pressure is at the extremely low level of about 0.2 psi, which 
is probably insufficient to support any appreciable heat transfer to components 
located within the compartment. 


20 


Compartment Volume = 3000 ft 



Skin Differential Pressure as a Function of Time for a Typical Descent Trajectory 



Compartment volume = 3000 ft 
4:1 elliptical orifice 



Fig. 11 - Compartment Pressure History During Typical Descent 


Compartment volume = 3000 ft 
4:1 elliptical orifice 
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Fig. 12- Compartment Temperature History During Typical Descent 


Section 5 

CONCLUSIONS AND RECOMMENDATIONS 

This investigation reveals, on the basis of computer calculated results, 
that shuttle compartment venting during the ascent phase of the trajectory is 
significantly more critical than compartment venting during the descent tra- 
jectory. This conclusion is based on a comparison of calculated compartment 
pressures for ascent and descent trajectories. The study further revealed 
that sufficient inflow venting is likely to occur during shuttle reentry on leak- 
age alone. This conclusion is based on the assumptions that the leakage rates 
of shuttle compartments are comparable to that on Saturn vehicles for a given 
compartment volume, and that a crush load differential pressure of 1.0 psi 
across compartment walls is within the spacecraft design stress limits. The 
vent ports required to vent the compartment gases during ascent, therefore, 
should be sufficient for any inflow venting requirements during descent. 

Although the computed compartment temperatures during descent do not 
appear to be excessive, it is recommended that the vent ports be located on 
the backside of the vehicle, as it is oriented during descent, where cooler ex- 
ternal air temperatures are expected. It is also recommended that further 
study be made of the severity of the thermal environment. 

An investigation should be made to determine the extent of structural 
vibration and panel flutter assuming the crush load pressure differential 
across the vehicle skin that is expected to occur during the descent trajec- 
tory. If the structural design requirements are such that a positive outward 
pressure differential across the vehicle skin must be maintained, then the 
compartments must be pressurized during reentry by some means. One 
method would be to pressurize the compartments by releasing gas from 
bottles stored on the shuttle vehicle. This method obviously results in a 
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severe payload penalty. Another method would be to utilize scoops for 
forcing air into the compartments. 

The compartment venting calculations presented in the example of 
this report are based on the assumption that the local external flow condi- 
tions are identical with the freestream conditions. Venting studies on the 
final shuttle design configuration should be based on actual local conditions 
at the vent port. These local conditions may be established by wind tunnel 
tests. 
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Appendix A 


INFLOW VENTING COMPUTER 
PROGRAM USER'S MANUAL 



Appendix A 


A. 1 APPLICATIONS 

The inflow venting program may be used for any situation of inflow 
venting, although it was specifically designed for descent trajectories of 
space vehicles. A single compartment with up to four external vents can be 
analyzed. Vents may have zero area, which can represent a location on the 
compartment wall where pressure differential calculations are needed. Also, 
one vent area change is allowed during the analysis. If the inflow discharge 
coefficient data input is not in a useful table form, a subroutine is called to 
interpolate and extrapolate the data into a uniform table. 

A. 2 INPUT REQUIREMENTS 

The input data consist of initial conditions, local parameter histories, 
and inflow orifice efficiency data. These are read in as follows: 


Card Set 


Variable 


Format 


1 NC, NV, NT 313 

where NC = 1 (single compartment) 

NV = number of vents (1,2,3 or 4) 

NT = number of local parameter time 
cards 

2 TT(L2), TP(L1,L2), 5E14.8 

TR(L1,L2), TV(L1,L2), 

TMACH(L1,L2) 

where TT(L2) = flight time, sec 

TP(L1, L2) = local static pressure, psia 
TR(L1, L2) = local density, slug/ft 3 
TV(L1, L2) = local velocity, ft/sec 
TMACH(L1,L2) = local Mach number 

NOTE: NT cards are read in for each of NV vents. 

All time points must be the same for each 
vent. 
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Card Set 


Variable 


Format 


3 


4 


5 


Ml, N, M 313 

where Ml = number of input table pressure 
ratios (and efficiencies) per 
Mach number 

N = number of Mach number inputs 
M = number of pressure ratios to be 
tabulated per Mach number 

NOTE: If M = 0, the input table is already of 
adequate length and in the required 
form with the same pressure ratios 
for each Mach number. Otherwise, 
a nonlinear interpolation develops a 
uniform pressure ratio table for all 
Mach numbers. 

DOX.DX2 2F10.5 

where DOX = initial pressure ratio of table 
DX2 = pressure ratio increment 

X1(I), PXZ(I,K),PY(1,K) 3F10.5 

where X 1 (I) = Mach number 

PXZ(I, K) = pressure ratio (compartment to external) 
PY(I, K) = orifice efficiency 


(5A) 
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XI (I), X2(NR), Y(I, NR) 3F10.5 

(replaces card sets 4 and 5 when M = 0) 

where XI (I) = Mach number 

X2(NR) = pressure ratio 
Y(I, NR) = orifice efficiency 

VOLL, PL, TL 3F10.5 

3 

where VOLL = compartment volume, ft 

PL = initial compartment pressure, psfa 
TL = initial compartment temperature, °R 
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AH(L1) 

2 

where AH(L1) = vent area, in 


E14.8 


» T, DT, TM, TC 4E14.8 

where T = starting time of analysis, sec 
DT = time increment, sec 
TM = stopping time, sec 
TC = time at which vent area change 
occurs, sec 

NOTE: T and TC must coincide with specific local data times. TM may 

be anytime during the flight. If no area change is used, TC may 
be left blank. 


» 
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A-2 



Card Set 


Variable 


Format 


9 KC.AH(KC) I2.E14.8 

(used only when area change is made) 

where KC = vent for which area is changed 
AH(KC) = new vent area, in^ 


A. 3 , RESTRICTIONS 

The program limitations basically involve the relation between compart- 
ment and external pressures. Analysis begins with compartment pressure 
either equal to external pressure (equilibrium), or less than external pressure; 
the latter case must be interpreted as the instantaneous opening of all vents 
at the starting time since the initial condition is zero venting everywhere. 

During the analysis, compartment pressure may not exceed any local 
external pressure. This is required since provision for outflow venting is 
not made in the program. Thus a short re -ascent in the trajectory is allow- 
able only if the pressure differential immediately prior to this re- ascent is 
large enough to ensure inflow throughout. This cannot be determined before- 
hand; if an outflow situation results from the trajectory, then the time point 
of changeover must be found and all data at this instant may be input as the 
initial conditions for the outflow program (Ref. 7). Future versions of the 
venting program will encompass both inflow and outflow. 

The possibility of outflow venting may also exist for multiple vents 
which have considerably different local flow conditions; it is conceivable that 
different vents may involve inflow and outflow simultaneously. Again, no 

allowance is made for this and the outflow program must be coordinated with 
the analysis. 

A. 4 INPUT LISTING SAMPLE 

A sample input data set is presented below for a single-vent, constant 
area analysis. Section A. 2 describes these in detail. 



1 1 *»4 

• 2 1 ^noonn+r»3 
.77OO000o+03 

. 23000000+03 

• 2*6?000O4-03 
.26000000+03 
•26440000+03 
•36000000+03 
•37OO0009+O3 
•380000OO+03 
•3QOO00Or4O3 
• 3997000O+03 
.40000000+03 
,4 1 ^0000*4*3 
.47090009+03 

• 4 300 OOO 04 O 3 
.44000000403 

• 4 POO 00 00403 
.49000000+03 

• •>0000000403 
.5 ! 000000403 
.P20O00OO403 
•*3O0O0Op+03 

•**000000403 

•*60p000o403 

•6 672000 04 03 
•67000000403 

• *7 l 3000^403 

• *ROOOOO A403 

• *r>oppn00403 

• 6P0O0O0<%4O3 
.78oo00Oo+03 

.70^00000403 

.6^000000+03 

• 1 1^74000+04 

7? ! * O 


•20800000-03 

• 13800000-03 
.69000000-04 
.69000000-04 
•69000000-04 
•69000000-04 
•1 1390000-01 

• 19790000-01 
.29720000-01 

• 37780000—0 1 
• 404 8000 o—O 1 
.40480000-0 1 
.40 1 40000-0 1 
•39370000-01 
.39030000-01 
.40000000-0 1 
•42910000-01 
.95340000-01 

• 129*0000400 
.17940000400 

• 19508000400 
•27179000400 
.484300^0400 
•8! O6OOOO4OO 
.94940000400 

• 1 00700nn+o l 

• 1037000040 1 

• 1291000040 1 

• 1 7*1 000040 1 
.40877*0040 1 
•*1 1 3000040 1 
•*7*40000401 
.*4*8000040 1 
•1011 0000407 


.5764611 1-07 
•3*1*4894-07 
• 279R329*— 07 
.14409167-07 
.145*3837-07 

• 15084294-07 
.19771614-0* 
.34653370-0* 
•53865592-05 
•69*96081-05 
.7*151473-05 
.75151473-05 
.74399137-05 
.729138*5-0* 
.71949386-05 
.741 1944 1-05 
.79954470-05 

• 18961460-04 
.26465029-04 
.36890938-04 
.52085200-04 
•77720000-04 

• 13300077-03 
.17399070-03 
.20446747-03 
.21623394-03 
•22259944-03 
.27738178-03 
.37604645-03 
.85198700-03 
.10087381-02 
•10316975-02 

• 10636164-07 

• 17*7440 1-02 


• 1 1 3794004-05 

• 1 1 7899004-0* 

• 1 12677004-05 

• 1 12422004-05 

• 1 12399004-05 

• 1 12384004-05 

• 102076004-05 
.04 4790 00 4-0 4 

• 57*160004-04 

• 807520004-04 

• 744420004-04 

• 738930004-04 

• 667740004-04 
.591310004-04 

• 5101 1 0004-04 
.447340004-04 

• 387640004-04 

• 233470004-04 

• 206050004-04 

• 1 79n<50004-04 

• 1 56490004-04 

• 136979004-04 

• 101640004-04 

• 837700004-03 

• 61*600004-03 

• 661*00004-03 

• 68 1 800004-03 

• 834300004-03 
*946*00004-03 
.94 1487004-0 3 

• 406800004-03 

• 416100004-03 

• 437800004-03 

• 334300004-03 


• 124630004-02 

• 127720004-02 

• 127470004-02 

• 127180004-02 

• 127150004-02 

• 127140004-02 
.943400004*0 1 

• 883000004-0 1 
.829800004-0 1 
.773100004-01 

• 714800004-01 

• 709900004-01 

• 6404 00004-0 1 

• 567100004-01 

• 4891 00004-0 1 
.474500004-01 

• 373100004-0 1 
.232100004-01 
.206300004-0 1 

• 181 00QO0+0 1 
.159080004-0 1 

• 139*80004-01 

• 1 04000004-0 1 
.860000004-00 
.636000004-00 
•68700000+00 

• 704000004-00 

• 862000004-00 

• 978 00 00 0+00 

• 968402004-00 

• 403^00004-00 

• 41 1000004-00 

• 431000004-00 

• 31 0000004-00 


C *AMpt_r FOO ONF *»4CH NUMPFR ONLY) 


1 • ■f'OOO 
1 • *9ono 

1 . *nnnn 
1 .*oooo 
1 . «99oo 

1 . *nooo 

] . «nnpo 
l .*9900 
1 .*0099 
1 .*0909 
I . mc*r\fsfs 

1 .*9999 
1 . *9099 
1 .*0000 
] .*9990 
1 ,40000 
1 ,40000 

1 # *ooon 

1 ,*0000 

1 ,«OO0O 

1 ,400^0 
1 *«oooo 


. oooor 

,04000 

* 1 OOOO 
« 1 *onn 
.2Apo>o 
• ?*onn 

* 10000 
. 14000 
,4 n OOO 
•4*000 
.^''OOO 
. 4400 ft 
,fOOOO 
.6*000 
• 70000 
.7*000 
.80000 
• 8*000 
*00000 
* 0*000 
.08000 

1 .00999 


,44404 
.*6194 
.5*498 
.*4707 
. *3008 
.*3467 
. *7666 
•*1741 
.*067? 
.49309 
. 47004 , 

.46*98 

.44097 

.43799 

.40007 

• 38404 
.3*008 

• 3080* 
. 7309* 

• 1469 * 

.o-rooo 

.09999 


■40 00. 
*0. 9 
4 30.9 


*•616 4*7 . 

noo.o 
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Appendix B 
PROGRAM LISTING 



c fhutti r percent i nvfnt T ng 

r 

C GFNFPAL SINGLE COMPARTMENT, FOUR VFNT INVENTING PROGRAM 

r 

r OUTVnNTiNG I R NOT Al-I OWABL<~ 

r 

r onF \/fnt apfa changf If allowfp 

r 

c initial condition is zero venting everywhere 


r 

r 

PR TNT OQO 
ono FORMAT <1 HI ) 
print 003 

otp FORMAT ( 20H OFF! NIT I ON OF OUTPUT SYMBOLS// 

14H DNC 1 2X45HI NSTANTANEOUS MASS FLOW RATE INTO C0MPARTMENT3X 1 OH ( SLU 
?r/<rr ) / 

34H DM J 1 2X4 1 H I NST ANT ANEOUS MASS FLOW RATE THROUGH VENT3X 1 OH ( SLUG/SE 
40/ 

F2H K i AXAOHDIRCHAPGo coffficifnt /) 

PRINT 033 

OTT FORMAT (6H MACH J1 0X20HVFNT JFT MAC.H NUMBER/ 
pTH Mr 1 0X1 OHM ASS IN COMPARTMFNT3X7H ( SLUGS ) / 

B3H Pr] 3X2OHC0MPARTMENT PRESSURE 3XBH ( LB/ I N? ) / 

7BH M j/ML 1 3X42H JET TO LOCAL MASS FLOW PER UNIT ARE A PATIO) 

PRINT op/, 

OTA FORM AT (/ 

JTH Pi. 1 3X31 HLOC AL PRFSSURF AT VENT LOC AT I ON3X8H ( LB/ 1 N? ) / 

3SH Pr 1 1 X23HCOMP ARTMFNT GAS DFNS I TY3X 1 OH ( SLUG/FT 3 ) / 

3H TC 1 3X27HCOMPAPTMENT GAS TEMPERATURE 3X6H ( DFGR ) / 

=;oh T T MF 11X11 HFL I GHT T I MF3XBH ( SFC ) / 

6 1 OH VENT AREA6X 8HPER VENT * 3X5H(IN2)/7H VOLUME, 

74H 1 2X1 BHCOMPARTMFNT V0LijME3X5H ( FT3 ) ) 

D IMFMS ION TT ( BO ) , TP< 4 «S0 ) , TR <4»S0 ) , TV < 4 » BO ) , TMACH ( 4 , BO ) 

D I MENS ION BP ( 4 ) , RR( 4 ) »BV<4 ) ,RH( 4 ) , VH( 4 ) , AMACH ( 4 ) 
n I M cnS I ON AH ( A ) , CH ( 4 ) , XV(4) 

D I MFNS I ON SRH ( 4 ) , SVH ( 4 ) ♦ SEH ( 4 ) , WH ( 4 ) 

D T MFNR ION Y I 30 ,30 ) , XI ( 3R ) , X2 ( 30 ) 

D IMFMS I ON PY( 3B* 3B ) »PX2 t 30,38) , X2X < 30 ) « Y 1 ( 30 ) 
n I MF me I ON BTTl ( 4 ) ,RTT2( 4 ) ,OMJ( A ) 

r PfAo IN and POINT TPAJPTTOPY data 

PRINT ooo 

PRINT 7 EC 

700 FORMAT! 1BH TRAJECTORY DATA//) 
read pr 4 , NC,NV, NT 
PC4 FORMAT! 3 13) 

DO 1?0 |.1 = I,NV 1 

PP I nit top , L 1 

tbs FORMAT! / c H VENT,I3»26H LOCAL EXTERNAL ^ CCND I T I ONS// 1 2XRHT I ME • SECAX , 
1 1 3HPPF SSURF , PS I A3X 1 4HDFNS I TY , SL/FT32X 1 0HVFL0C I TY , FT /SF C6X 1 1 HMACH N 
•PI IMRFP/ ) 
no 1 00 L2 = 1 , NT 

READ 80 1 , TT ! L2 ) « TP ! L 1 ,LP ) • TR !L1 »L2) , TV < L 1 , L2 ) , TMACH (LI , L2 ) 
po 1 format ( or i 4 ,p ) 

WRITF(6«760) TT(L2> ,TP(L1 ,L2> , TR ( L 1 ,L2) ,TV(L1 ,L2) , TMACH ( L 1 , L2 ) 

TFO PORMAT ( pyoF 1 7 , b ) 

120 TP(L 1 «L? ) =TP(L 1 ,L2)* 1 A4 • 

c read in and print venting efficiency data 
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PF An RFA.M1 ,M,M 
!F(M) 

P 7 A M=M| 

no flpp 1 = 1 « M 
no n-sn mp= i «m 

arr prAn rrr.xi ( n .xrinri »v< i »nr> 

FORMAT ( FFl n ,c ) 

OO Tn flOl 
m? CONT INI if 

pfao Rnc,nnx,nx? 
no Ann I = i , n 

00 7 ic: 1 ,M 1 

RFAn flf'F.vi ( n ,py?t I ,<) ,PY(! ,K) 

VI (K)rPV( f ,«) 

T VPV(U') =Dy? ( I » K ) 

PTAPTrOPy 
no i n mP= l , m 

y?(MO) s RTAPT 

CALL TAR I NT ( 1 « 7 » 1 .Ml , X?X,XP(NR) ) 

CALL TARINT ( 7. 7, 0. 0 . Y 1 « Y( I ,NR) > 
in STAp T =FTART+nXP 

AOn CONT TNI IF 
A O 1 CONTTMtlF 
op|mt non 
PRINT OPft 

P2F. FORMAT ( F 1 H TABLF OF DISCHARGE COEFFICIENTS 

PRINT RRP 

RRR FORMAT) 1 H ) 
no R I s 1 , M 
PPTNt opn « X 1 ( T ) 

F6RMAT(//6H MACHL2XF 1 O. 3 « 8XFHPC /PLQX4HK /) 

no R NP= 1 « M 

R PRINT on 1 , XP ( NP ) « Y ( I ,NP ) 

02 1 FORMAT ( 3?XF 1 O . 7, 3XF1 P .F ) 

OR tmt non 
0= 1 .A 

PRNr ( I . + ( G -1 • ) /?* >** (-G / ( O -1.)) 

OC= T T1 F. 

FYP= ( 1 .-0/0 
p PF An ( c , roc ) VOLL.PL.TL 
print o 1 n.vni.L 

° 1 n FORMAT) 1 OH COMPARTMENT VOLUMFFXF 1 0 •*=:/ ) 

PL = P| /(GC *TL ) 

WL=R| *VOLI. 
no 1 nc L 1 = 1 ,NV 
PFAO Rn 1 « A M ( L_ 1 ) 

PRINT OINLl.AH(U) 

0.7 1 FORMAT ( lOH VFNT ARFAI?«7XE10*F) 

IOc, AH(Li ) = AM(L1 )/1 Aa. 

PRINT RRR 
ORIMt RRP 

C RFAn in TIMF INFORMATION 

pfao not ,T«OT,TM,Tr 
no R 07 L. K = 1 .NT 

IF(T-TT)LR) )RO7,R0R,RO7 
P 07 CONTINUF 
RRR CONT TNI if 

C CREATF AND PRINT INITIAL OUTPUT 
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nniKjr nric . t 

OOB FORMAT < /'■-•H T1MFV4YF1 2.B) 

PL I - PI / 1 44 . 

OWt = , n 

1 0 PP I NT OOis, DL I . PL * TL « WL « OVL 

O0ft FORMAT ( 3H PC.6XEl?.B.3X2HRC6XE12.5.3X2HTC6XEl?.B,3X2HMC6XE12.r> 
1 . ByBMOMCByFl 2 ,b> 
print bpo 

STMi -r, 
flHT=A, 

no 111 L. 1 s 1 . MV 
PP(Li ) = TP ( L 1 (I.F) 

BR(L1 )=TR(L1 .LB) 

PV(I_1 )=T\/(L1 .LB) 

P T T 1 (L 1 ) =BP ( L 1 ) /GC/ER (LI )*( 1 •+( G“ 1 • > X 2 • *TMACH < L 1 * LB )**?.) 

BTM 1=BTM1+RTT1 (Ll )*AH(L1) 

AHT= AHT+AH(L 1 ) 

PH(|_ 1 )=PL 
xvd.n = ,n 
vhil 1 ) = .n 
P = RP(I 1 ) /I 44. 

PP= ( n° ( L 1 ) -PL ) XI 44 . 

IF< AH(L 1 ) > 1 00 . 1 00 . 1 0 1 
inn print ooi.li.p.DP 

R0 1 FORMAT (RX5H VFNT I 2 , 8X2HPL6XE 1 2 . B , 3X5HPL-PC. 3XE 1 2 . 5 ) 

GO T O 111 
1 O 1 D = PL /RP(L 1 ) 

CALL P! VLP( IP. TMACH ( L 1 . LB ) « X I . N ,R « X2 ,M . V * X< ) 

F H(u )=y« 

print pop.li .p.dp.p 

R02 FORMAT <9X5H VENT l 2 » 8X2HPL6XE 12.5* 3X5HPL— PC3XE 12.5* 3X5HPC/PL3X 
1F12.R) 

BIO RM=,n 

XL=TMAOH(L 1 (LR ) 

PRINT on b * XL * RM . y< 

903 FORMAT ( 24X5HM ACHL3XE 12.5. 3X5HM J/ML3XE 1 2 • 5 » 3X 

1 1 HK7yF1 2.B) 

1 1 n print on a ,yv(L 1 ) «OU'L 

904 FORMAT < 24XBHMACHJ3XF12.B.3X3HDMJ5XE12.B) 

1 1 1 PONT TNI IP 

I F ( Amt ) 31 1 , B 1 t . B IP 

31 1 TTM 1 -O, 
r.n ta b 1 

312 TTM 1 =FTM 1 /AHT 
BIB CONTI M( IF 

INTERPOLATE BOUNDARY CONDITIONS FOR NEXT TIME POINT 
J = 3 

k = o 

1 TrT+AT 
FS=K+1 

OO lf5 11 = KS . NT 

I F ( T_TT < L 1 ))14.112,11? 

1 1 ? K=K + 1 
113 CONJTTMUP 
14 FONT TNI IF 

IF(K-< J+l )/2) 16. IB. IB 

IP TF( K4-J/2-NT ) 1 7. 1 P . 1 P, 

1 6 N7 = 1 


l y 

1 * 
1 n 


? 1 4 


P o O 

PO 1 


pop 

1 1 c 


1 1 * 
] 1 7 

7 1 4 

71 « 
71 * 


O a 


4 1 


7« 


“7 "7 


no TO ]R 
N7=K-( J-l ) /? 
no TO to 
M7=MT-J+ 1 
FONT T Ml IF 

HO 1 1 4 t 1 sr 1 « MV 

amatm ( i i > = # r 

non t ) = 

HQU 1 > = .0 
PV(L 1 ) = .0 
PRINT RRS 
no n R I 7« 1 ♦ J 
Ml 1 SM7 + L1- 1 
0=1 . 

no 701 L4»1«J 
ML?=m7+L4- 1 

T F ( M| !-NL7)70O#?Ol *700 

C-C- ( T-TT ( ML 7 / 5 / ( TT ( ML 1 ) -TTC NL P) ) 

continue 

no PAP L 1 = 1 f MV 

AMA0H(L1 )=AMACH(L1 )+C*TMACH(Ll tNLl ) 

RP ( L 1 )=RP(L1 )+C*TP(Ll *NLl ) 

RR(L1 ) = RR ( L 1 > + C*TR ( L 1 % ML t ) 

PV(Li )=RV(Ll )+C*TV<Ll *NL1 > 

COMT TMUF 
CONTTMt IF 
o = i oa # 
e;TMP-n # 

no 1|7 L1=1|MV 

RTTPd 1 ) =RP ( L 1 )/rO/RR(Ll )*< 1 •+<G-1 • ) /2 • * AM ACH ( L 1 >**7) 
CJM9-CTMP+PTT P ( L 1 )*AH(L1 ) 

I F ( P— RP (LI ) ) 1 17# 11 7* 1 16 

P = RP ( I 1 ) 

FONT T Ml IF 

I F ( AMT ) T1 4 t 7 1 4 t 7 1 R 
TTMpr n # 

no To P1A 
TTM? = cxm?/AHT 
FONT T Mi IF 

• SFT iiP PQFSSUPF O IFFEPFNT I ALS AMD START TNG PRESSURE 
nPOr ( n-PL ) /7. 
np(_ t -nnn 
PL T =nl 
k 1 = 1 

FOMTTK'MF 

nwi •= # o 

■iterate for compartment venting solution 

nvi_ = nu/i i 

no i (i = i ,wv 

T F ( Am ( L 1 ) ) p7, 1 7P,7R 

PHI r=Ol. T 

d = pi T /non i i 

IF ( P-1 . )7ft*7(S«77 

oh T - pn ( | n V 

P= 1 . 

FOMTTMllF ' 

CALL RIVLP( I R » AM A CH ( L 1 ) ♦ X 1 ♦ N % R « X2 t M , Y • XK ) 

IF(P.PPN) 1 1R» 11 <p* 1 1*P 


V 


r 


4 
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1 1 - H I = HP ( L 1 “PPM 

1 1 o PRrPwt /HP (LI ) 

DRp=np**FXP 

PHI SDH I /r-C /HTTP (L 1 ) *ppp 

XV ( L 1 ) = SORT ( 2 • / ( G- 1 • ) * ( PRR— 1 . ) ) 

vh!=vv(li >*sopt(g*oc*«3tt2(li )/prr) 

CRH(| 1 ) = RHT 
5 VH ( | 1 ) = VH T 
PHH(( 11=XK 

HMJ(, 1 )=PH1*VHI«XX*AH(L1 ) 

WH( L 1 ) = < (RH(L 1 ) *VH( L 1 >*FH( L 1 ) ) +( RH I *VHI *XK 1 )*AH(L1 ) *DT /?. 
PW|. =p»*<I +WH(L 1 ) 

COMTTNHIF 

PPL = PL I - (PL+OWL*G*GC/2./VOLL* ( TTM 1 +TTM2 ) ) 

GO TP «7) ,«] 

AS !P( Apf ( PPL ) — . OO 1 44 ) 47* 47« 40 
A ~y GO Tp (F7,4P],i<| 

AO K!=0 

go tp h? 

ao go to ( ho « r i ) , < i 
ho <I=? 

PL.1C-DL1 
OPLP=PPL 
PL I =oi I +PPL T 
ro to on 

51 IF (APS (OPLM/DPLS+ABS ( DPL ) /PPL ) 5.3 . 52 . 53 
5? PPL I = PPL* ( PL I -PL I S ) /( DPLS-OPL ) 

PL I S = OL I 
pL.i=oi.n-npi i 
PPL T-npi. T/2. 

PPL C - = PPL 
G O To on 

52 I F ( AoS (OPLR)-ABS ( PPL ) ) 54 « 57 « 50 
ha PPL. I=-pPLI/2. 

GO TO =P 

5H IF(APS (PPLF)-ARS (PPL) )5B«57«57 
HS <1=4 

PL I =ni I s 
HP To P0 

5FT OATA fop MFXT I MTHPVAL 

H 7 PONT I Ml IF 
PL=P( I 
1«'L = W|+PWL 

TTM 1 -TTM? 

PMp=o . 

PO 1 -ja 1.1 = 1, MV 
IF(Aw(L1)> 1 .24 « 1 34 * 1 .20 
1 30 PH (Li ) = SRH (Li ) 

VH (Li ) = SVH (LI ) 

FH(L 1 ) = SFH (LI ) 
wrmMr+nMjiLi 1 
t HA C OMT T Ml l<— 

PRINT Ol IT d UT PAT a 

hit hRimt oos , T 

(SH R = liH /VOLL 

TL=P| /(GO *R> 

RI=P| / 1 44 • 

OR I NT oo* , D I , P , Tl. • WL « owe 
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PRIMt flRfl 

no i /i o l 1 = J «mv 
P = BP ( | i )/UR. 
np=(pp(| j )_dl )/i«4, 

TP ( Am ( I 1)) 107,lo7,tBfl 
1 Vr PPImt oo t , L) ,p , pp 
OO Tr> 1 A 0 

t to roMt r m( ir 

p = °l /nril I ) 

PR I n't nn?,| i ,n,np,p 

1 RM = PH( L1 ) *VH(L i ) /<RR(L1 > *PV(L1 > >*FH(Ll) 

P° I Nt om, AMAPH(L1 ) ,RM,FFH(L1 ) 

P c PRINT 0O4 « XV (Li ) , om J (LI) 

140 PONT T mi if 

IF ( T-TP ) 141)1 42 ,14 1 

140 PPAp O^B«FP «AH(Kp ) 

°0R FOPVATI l?)F14.fl) 

pqtmt nno.k'r . AWfi^r ^ 

°°° FORMAT ( ///RX1 FH. ... .VENT AREA I2»12H CHANGED TO E12.F,6H SQ 
AHfk-rirAMKPl/Hd, 

AHT=P . 

no 1 AO (1 = 1 ,MV 
140 AHT= AMT+flH f L 1 ) 

14 1 PONT T Ml ic 

IF(T-TM) 10,BF,AF 
AN P ONT I Ml IP 
PPImt non 
non cjop 
rNtn 


IN//) 


V 


r 


i 
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r 

c 

c 

c 

r 

c 

c 

r 

r 

r 

r 

v 


* 


* r 


♦ 


SUBROUTINE BIVLP( I FRR »X « XT ,N « Y, yT . M , 7T • 7 ) 

this SUBROUTINE uses table LOOKUP to find a dependent 

VARIABLE as a function of two independent variables. 

I FpR is an FRROR FLAG, WHEN NOT FQuAL 7 ERO = ERROR 

Y IS ARGUMENT OF FIRST INDFPFNDFNT VARIABLE 

YT IF MATRIX OF FIRST INDEPENDENT VARIABLE X ( I > • I 

N IS NO. VALUES OF FIRST INDEPENDENT VARIABLE 

Y P ARGUMFNT OF SECOND INDFPFNDFNT VAR. 

YT IS MATRIX OF SFCOND INDEPENDENT VAR. Y(J),J= 

M IS NO. VALUFS OF SFCOND INDFPFNDFNT VAR. 

7T IS MATRIX OF DEPFNDFNT VAR. 7<I«J>=7< X(I). Y(J) 

7 IS DFPFNOFNT VAR. ANSWER 

0 T MFms T ON XT ( SF) , YT(3S) , 7T(3B,SB) 

IF( VT( 1 )-X )3,S,pE 

o CONTINUF 

IF(YT< 1 )-Y)4«4»?E 
A CONT T Ml |F 
DO B I = P , M 

IF ( XT ( I )-X)B. 10,10 
P CONTINUF 

I F ( XT ( N ) -X ) PB , 101 , 101 
101 t=N 
1 O NH I = T 
Nt 0= T - 1 
Y 1 =XT (MLO) 
y?sYT(MHI ) 

no p b j=?,m 

1 F ( YT ( J > — Y ) PP. .16.1ft 
?P CONTTMt F 

I F ( yt ( m ) -v ) ps, j on, i rn 
1 O J = M 
1 Ei MHT=u 
MLO= J - 1 
yi =yt < v|. r , ) 

Y? = YT (MH1 ) 

70flr7T( NLO , MLO ) 

701 r7T ( MLO , MH T ) 

71 n=7T(NHl , MLO ) 

71 1 =7T(MMt ,MHI ) 

V1=7aO + ( 700-7 t O ) * ( Y-y 1 ) / ( Yl -YP ) 

\t? = 701 +(701-71 1 )*(X-X1 ) / ( X 1 —X? ) 

7 =V| 4-(\'1-V?)*( Y-Y1 )/( Yl-Y?) 

I rpp-n 
jom 

PS I FPP = } 

PQPOD f m DATA 

PFTi inK.i 
r no 


= 1 f N 
1 

) 
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C 
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c 

r 

C 

r 

r 

c 

r 

c 

r 

C 

r 

c 

c 

c 

c 

r 

r-- 


r ali 

KK 


NTC *^=1,4,5 1 OP GREATFR. 

0 OR LFSS* 


NTE 

TAB 


KK= 1,4,5 
R'Ks.l* 

KK= 1,4,5 

<K = P 


P 

o 

QOn 


1 n 

1 On 

1 1 

1 n 


SUBROUTINE TARINT ( KK ,N I P , NTS , NTE , TAB , ARG ) 

DETERMINES from tabled data the correlative VALUE oe 

the prescribed ARGUMENT by means of the lagrangf inter- 
polation equation 

TAB I NT (KK « NlP»NTS,NTF , TAB , ARG) 

K'ONTROL konst ant « 1 — SEARCH* COMPUTES COEFFICIENTS 

3“ APPLY, COMPUTES INTERPOLATED VALUES 

s-fxtract, outputs coefficients 

4~ SUPPLY • COMPUTES COEFFICIENTS FOR MISSING 
TARLF VALUES 

S-7ER0* COMPUTES COEFFICIENTS FOR ZERO VALUE 
NIP KK = - S NUMRER interpolation points to be useo. 2 OR GREATER 

BUT LFSS THAN OR EOUAL TO to OR (NTF-NTS+1) WHICHEVER 
IS SMALLER 

INITIAL TARLE INDEX AND STARTING POINT 
OF SEARCH 

(MUST FOLLOW AN INITIAL CALL WITH NTS=1 
OR GRF A TER • ) iNTIAL TAHLt: POSITION RE- 
MAINS AS PREVIOUSLY DEFINED RUT SEARCH 
RESUMES AT TAELS' INDEX OE LAST SEARCH 
TERMINATION 

TERMINAL TABLE INDEX AND CUT-OFF POINT OF SEARCH 
TARLE INDEX MINUS ONF 

INDEPENDENT TABLE VARIABLE (SINGLY INDEXED FOR TABLE 
POSITION) 

DEPENDENT table VARIABLE (SINGLY INDEXFD FOR TABLE 
POS I T I ON ) 

NAME DESIRED FOR COEFFICIENTS (SINGLY INDEXED FOR 
COEFFICIENT NUMBER) 

ARGUMFNT OF INDEPENDENT VARIARLF 

NAME DFSIRRD FOR INTERPOLATED DFPFNDFNT VARIARLF 
VALUE OE INDEPENDENT TABLE VARIABLE FOR WHICH DEPEN- 
DENT TABLE VARIARLF IS TO Br COMPUTED 
OUTPUT DATA 

TAB ( 1 ) , COFF (10) 

GO TO ( 1 . !7,?n,i ,??) ,«« 

-ffARoh and SUPPLY PORTION of interpolation 

IP (N TP. GT. NTF-NTS+1 ) GO TO 99 
c IGM-1 . 

IF (N|T«; ) /, % & , 

NT SS-NTS 

( 1?:MTC 

GO TO R 
MTF = | 1 R 

TF(TAR(MTSS)-TAB (NTSS+1 ) ) 7,7,* 

F I GM — — ) . 

OTAr = . 1 * R I C-N* ( TAR ( NTF ) -TAB ( NTSF ) ) 

IF { R TON* < APO- ( TAR ( NTFR ) -OTAR ) ) ) O, R , p 
IF ( F f AN* ( ARG- ( TAB ( NTF ) +DT AR ) ) ) 10,10,9 

">P I Te* ( a , nrn ) ARA , T AB ( NTSS ) ,TAR(NTF) 

I b™ , ™;m'?^p H c^' DE " aNGE ° F "BUP.EI6.8.3H T0.E16.8.25H 
no 1 no |_ j =t ] R , fv T p 

I F ( r | GN* ( ARG-T AR (L 1 ) ) ) 1 1 , t OO , 1 00 

CONT T Ml IF 
L 1 F = t 1 - 1 

IE(L 1 F-(NIP+1 ) /? ) 13,12,12 

TF(LlF+NTP/p-NTF) 14 ,ir,ir 


ARG 


KK=3* 

KK = 1 
K 'K = P* 
K K = 4 * 


-+ DFmOxpr 
O I MFMS I OKJ 
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m 


* 


4 


m TO 

14 NSP = | 1 «'.-(NIP-l )/P 

IF ( mco— NT C '^ ) 10 . 16 , 1 * 

IF NFP=mTF-NIP+1 
IF. HO Ini L 1 = 1 « N I P 
ML. 1 =mfP+L 1 -1 
0=1 . 

no PAO L?= 1 ,MIP 

NL2=NSP+L?- 1 
IF(|_i •ro.L?)r.O TO pop 
I F ( ApO-T AR ( NLF ) ) POp , pn j , pop 
00 1 TF (K , k',=0«4)AO TO ?n.T 

POP 0 = C* ( AOO-f AP(ML?) )/( TAP (ML 1 ) -TABLNL? ) ) 
pop CONTTNUF 
in i ro<=F(i i )=r 

k'Cs 1 
RETURN 

APPLv ppp-fiON OF INTERPOLATION 
1 I F < <n,AT *0)00 TO IP 

IP WP t Tc ( f, « POP ) 

BOO FORMAT ( 1 H ,On H COFFFICIENTS WERE NOT COMPUTED) 


IP ARAr.n 

DO ]n? L. 1 = 1 « N I P 
ML 1 =mcp+l 1 - 1 

I"'? ARO= AQO+OOFF ( L 1 ) *TAR(NL1 ) 

PRTl inM 

FVTRACT intfrpolation COFFFICIFNTF 
Po IF(K'r) 1R* 1 R ,?1 
P\ no |m L I = 1 »NI d 
lOo TAR (| 1)=C0 Ff<L1) 

NTF = m<;p- ] 

P^TI |on 

7FPO TAPLP nOFITION 
P P APO,= .P 

nep=mtc- i 

no 1 PP 11=1 ,NFP 

IF(TAR(L1)/APF<TAR(L1 ) )-TAP(Ll+l )/ABS(TAR(Ll+ 
in* LI F = | 1 +1 
no to i i 
IOr CONTINUE 

WR I Tc - ( F« pnj > 

FOi FORMAT! 1H ,??HNO 7FR0 TARLF POSITION) 

kro = o 


) ) ) 10F, 1 04 , I 05 


PFTl 'DM 

°° WP I Tr~ ( A * ooo ) 

ORQ FORM AT ( 1 HI « 68HNUMBER OF 
1 FF TNTFRPOLATFD ) 
cron 

PNO 


INTERPOLATION POINTS GREATER THAN TABLE TO 
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